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Abstract. Consider a low temperature stochastic Ising model in the phase coexistence 
regime with Markov semigroup Pt. A fundamental and still largely open problem is the 
understanding of the long time behavior of 5,,Pt when the initial configuration r) is sampled 
from a highly disordered state v (e.g. a product Bernoulli measure or a high temperature 
Gibbs measure). Exploiting recent progresses in the analysis of the mixing time of Monte 
Carlo Markov chains for discrete spin models on a regular b-ary tree T*", we tackle the 
above problem for the Ising and hard core gas (independent sets) models on T*. If 
is a biased product Bernoulli law then, under various assumptions on the bias and on 
the thermodjTiamic parameters, we prove //-almost sure weak convergence of 5r,Pt to an 
extremal Gibbs measure (pure phase) and show that the limit is approached at least as 
fast as a stretched exponential of the time t. In the context of randomized algorithms and 
if one considers the Glauber dynamics on a large, finite tree, our results prove fast local 
relaxation to equilibrium on time scales much smaller than the true mixing time, provided 
that the starting point of the chain is not taken as the worst one but it is rather sampled 
from a suitable distribution. 



1. Introduction 

Let G = {y, E) be a countable infinite graph of bounded degree and consider, for def- 
initeness, a continuous time stochastic Ising model (Glauber dynamics) {o'I'}(>g on G 
with initial condition rj and infinitesimal generator C. Here rj is picked from the set of 
assignments of a ±1 variable to each vertex x & V . The main problems discussed in this 
paper can be formulated as follows. 

Assume that the thermodynamic parameters are such that there exist multiple reversible 
Gibbs measures for C For stochastic Ising models this amounts to say that the inverse 
temperature /? and the external field h are such that //+ ^ fi^ , where /i+ and //~ are 
the Gibbs measures obtained by taking infinite volume limits with pure + and — boundary 
conditions, respectively. Suppose that G is distributed according to a Bernoulli product 
measure with parameter p, i.e. {rix}x&v is a collection of i.i.d. random variables with 
P(ry^ = +1) =p. Then: 

i) Under which condition on the bias p is the Ising plus phase ^+ the unique limit point 
of the law of a^ast^ oo, for a.a. rj ? 

ii) If so, how fast does the law of approach ^+ ? 

The above questions, with G some regular lattice, have their origin in the theory of 
"phase ordering kinetics" [jSj - that is growth of order through a dynamical domain coars- 
ening - and clearly represent basic problems in the theory of interacting particle systems. 
Unfortunately, a rigorous approach these problems is still largely missing. 

If the law of the starting configuration r] stochastically dominates the plus phase /i+ it 
is possible to use some monotonicity arguments (allowed by the ferromagnetic character 
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of the model) to prove that /i+ is indeed the unique hmiting point of the process and that 
the convergence takes place faster than any inverse power of t. We refer to section 6.7 in 
||T8 ] for the case G = Z'^ and to our Lemma lT4l beIow for a stronger statement in the case 
of regular trees. However, it is easily seen that such a stochastic domination requirement 
for the initial Bernoulli distribution forces the bias p to be exponentially close to 1 when 
/? ^ oo. For (3 < oo we do not know of any result that goes beyond this simple case. 

On the other hand, the extreme case /3 = oo (zero temperature Glauber dynamics) has 
received considerable attention in the probabilistic literature, and various kinds of graphs 
(Z'^, the hexagonal lattice and the binary tree) have been considered [28, 30, 29. 31[l8ll. 
In this case, besides the motivation from physics to study simple models of spatial domain 
coarsening, there is also an interesting connection with (non-linear) voter models 111 711 . 
The relevant quantities are then the probability that a given vertex flips its value finitely 
or infinitely many times, the probability that a given spin has not flipped before time t, the 
typical size of clusters of vertices with a common spin value and other related percolation 
questions. 

Going back to our original problems, a major obstacle for progresses in the case = 1/^ 
is represented by the absence of tight bounds on the mixing time of the Glauber dynamics 
in finite boxes with plus boundary conditions, i.e. those boundary conditions that select 
the plus phase. On the contrary, when G is the regular 6-ary tree, this question has been 
recently solved in a sharp and constructive way for various models [19]. Exploiting the 
results of [.191 we have been able to study the above basic questions for two attractive 
systems on trees: the Ising model and the hard core gas (independent sets). Our results 
provide some answers to i) and ii) in non trivial cases. For instance we show that if the 
bias is sufficiently large (but independent of /3) then we have the desired convergence for 
all temperatures. The paper also includes a discussion of several interesting problems that 
are left unsolved and that we would like to consider in future work. For simplicity we 
present now our main result only for the Ising model, and defer the reader to section 6 
for the analogous theorem for the hard-core gas. Before stating our results we will now 
briefly overview the model and its basic features. 

1.1. The Ising model on the 6-ary tree. From now on denotes the infinite, rooted 
6-ary tree, where each vertex has exactly h children (6 ^ 2 is a given integer) . The Ising 
Gibbs measure on at inverse temperature {3 and external field h, formally given by 



where E is the set of edges of T'', has recently received a lot of attention as the canon- 
ical example of a statistical physics model on a "non-amenable" graph (i.e., one whose 
boundary is of comparable size to its volume) - see e.g. Il3l 1121 171 1261 II l[ in|2l l . The phase 
diagram of the model in the (/i, /3) plane is known ([l9l ll4ll ) to be quite different from that 
on the cubic lattice 1/^ (see Fig. ^ . 

We now recall some of its basic features. We write Ti for the rooted tree obtained by 
removing all vertices which are at distance greater than I from the root. The measures ^+ 
and are obtained by imposing + 1 and, respectively, — 1 boundary data at the leaves of 
Ti and taking the limit I ^ oo. The free measure is defined as the limit I ^ co when 
the boundary data at the leaves of are free (i.e. absent). 
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Figure 1. The critical field hc{/3,b). The Gibbs measure is unique above 
the curve. 

On the line h = there is a first critical value /3o = ^log[j^), marking the dividing 
line between uniqueness and non-uniqueness of the Gibbs measure (i.e. /i+ / as soon 
as (3 ^ Pq). Then, in sharp contrast to the model on Z'^, there is a second critical point 
I3i = \ log(^ii) which is often referred to as the "spin-glass critical point" (SI and has 
different interpretations. If one considers for instance the model with /i = on the finite 
tree with i.i.d. Bernoulli random boundary data r/ with p = 1/2 at the leaves of T^, 
then the distribution of the magnetization at the root (as a function of rj) becomes non 
trivial only if /? > see [6]. In particular, as £ ^ oo, for j3 ^ Pi the Gibbs measure 
on T£ with the above random boundary -q, converges (weakly) a.s. to the free measure 
^free Another way to look at /3i is to say that j.i^^^'^ is an extremal Gibbs measure iff 
P ^ Pi (see [3I II21IT3I ITII and, more recently, [19]). Finally Pi has also the interpretation 
of the non-reconstruction/reconstruction threshold in the context of "bit reconstruction 
problems" on a noisy symmetric channel (71 1221 12 111 . 

When an external field h is added to the system, it turns out that for all P > Pq, there 
is a critical value he = hdP, 6) > of the field such that /i+ 7^ /^^ iff \h\ ^ he. The Ising 
model on the tree at external field h = ±hc therefore shares the following two properties 
with the classical Ising model on at zero external field: on one hand the Gibbs measure 
is sensitive to the choice of boundary condition; on the other hand any arbitrarily small 
increase of \h\ causes the Gibbs measure to become insensitive to the boundary condition. 

1.2. The Glauber dynamics. The Glauber dynamics on T'' is the unique Markov process 
Wt}t ^ on with (tJLq = rj and Markov generator C formally given by 



where denotes the configuration obtained from a by flipping the spin at x, and Cx{cr) 
denotes the flip rate at x. 

Glauber dynamics on trees has received recently considerable interest 0] I19 II . Results 
in [19] show in a rather strong form that the mixing time (see e.g. II24II for a definition) 
on the finite subtree is always 0{i) if either P < Pi and h is arbitrary or if P,h are 
arbitrary and the boundary conditions on the leaves of Te are identically equal to -l-l (or, 
by symmetry, to —1). In particular, the Glauber dynamics in the pure plus phase ^+ always 
mixes fast (see II 1911 and section 2 below for more details) . 




(1.1) 
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Although all our results apply to any choice of finite-range, uniformly positive, bounded 
and attractive flip rates satisfying the detailed balance condition w.nt. the Ising Gibbs 
measure (see [ 18 ]), for simplicity in the sequel we will work with a specific choice known 
as the heat-bath dynamics (see section 2 below for the definition) . We will use the standard 
notation Pt = e*^ for the Markov semigroup associated to £. The spin at x at time t with 
starting configuration i] is denoted by a^^ and we will often use the shortcut notation 

Pt,-Av) = {Pta^M (1-2) 

for the expected value E((tJ'^) of a^^ given that the process starts in rj. 

1.3. Main results. In order to state our main results we need an extra bit of notation. We 
first define the set of initial configurations rj such that the Glauber dynamics converges 
weakly, at a certain rate, to the plus phase 

Definition 1.1. Given a G (0, 1), S7q will denote the set of starting configurations G $7 such 
that for any x G there exists a time to = to(??, x) < oo such that for all t ^ to 

I PtAv) - l^^i^x) I < exp(-t") . (1.3) 

We will see in Corollary 12.31 below that for any ry G f^a the law of the process 
converges weakly to /i+ as t ^ cx). 

The initial configuration rj is often sampled from a Bernoulli product measure with 
parameter p, i.e. rj^ = +1 with probability p and rj^ = —1 with probability 1 — p indepen- 
dently for each x G T*. We write Pp, Ej, for the corresponding probability and expectation. 
Finally, we need to recall the notion of partial ordering (stochastic domination) between 
probability measures on Q. Given two configurations a,r] G $7 we will write u ^ r/ iff 
a{x) ^ r/(x) Vx G T''. A function / : M is called monotone increasing (decreasing) if 
a a' implies f{a) ^ /(o"') (/(o") > /(o"'))- Given two probability measures fi, fi' on Cl 
we write /i < fi' if fi{f) ^ ^'(/) for all (bounded and measurable) increasing functions /. 

Our main results can now be stated as follows. 
Theorem 1.2. 

a) For every a > 0, b ^ 2, there exists p < I such that for all /? G (0, oo) and h ^ —hdP, b)+a 
we have ^{Qa) = h for some a = a(/3, h, b) > 0, for any initial distribution v such that 

b) For every p > \, there exist 6o £ N and Po G (0, oo) such that for h = 0, b ^ bo, P ^ /3o 
we have ^{fla) = I, for some a = a{P, b,p) > 0, for any initial distribution v such that 

c) Let denote the event defined in 1L3P with in place of /x"*". For every p < 1 there 
exist 5o G N and Pq g (0, oo) such that for b ^ bo, P ^ Po and h = —hc{P,b), we have 

= l/o'" some a = a{P, b,p) > 0, for all initial distributions v such that v ^ Pp. 

1.4. Remarks. Let us make some remarks on the above statements. 

1. We believe that in the case /i = 0, convergence to the plus phase should occur as soon as 
p> \- Unfortunately our bounds on p in statement a) are far from being sharp. However, 
as stated in bj, we can approach the critical value \, by taking b large. Another interesting 
issue is the dependence of p on h. Our technique in the proof of part a) of Theorem 11.21 
breaks down in the case h = —hdP, b) and the value of p in that statement approaches 1 
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as a ^ 0. On the other hand, statement c) shows that if /i = —hc{P, b) the critical value of 
p for convergence to the plus phase must approach 1 when 6 ^ oo. 

2. The main arguments we use to prove Theorem ll.2l are based on two essential features of 
the Ising model on T^. The first is monotonicity which is shared by all so-called attractive 
interacting particle systems. The second is the so-called "rigidity" of critical phases for 
spin systems on trees ([2]). Roughly speaking the latter means that, as long as h > —he, 
if we are in the pure phase ^+ we can add a small density of spins of the opposite (— ) 
phase and this will not alter significantly the structure of This, in a sense, is what we 
do when we introduce obstacles (see section 3 below) to lower bound the magnetization 
Pt,x{'ri)- The hard core gas model will be shown to have both these properties and our 
results there (see Theorem 16.11 below) will be obtained essentially by the same methods. 
On the other hand these techniques do not apply when there is no rigidity of phases, as 
e.g. in the Ising model on Z*^ (see II25II for a deep investigation of the metastable behavior 
of this model), or when there is no attractivity, as e.g. in the g-state Potts model for (7^3. 

3. A close check of the various probabilistic estimates needed for the proof of Theorem 
II. 21 and which are behind a Borel-Cantelli characterization of the set JIq, shows that there 
is also an -version of Theorem ll.21 with (II. 3D replaced by a bound of the form: 



for any t large enough. 

1.5. Plan of the paper. The rest of the paper is organized as follows. In section 2 we 
give the basic preliminaries for the proof of Theorem 11.21 In section 3 we explain our 
main argument. In particular, here we give the proof of Theorem 11.21 by assuming the 
validity of several technical claims. Section 4 and 5 deal with the proof of these claims. 
In section 6 we present our results for the hard core gas. Finally, some further results and 
open problems will be discussed in section 7. 



Here we first collect several useful preliminaries concerning the Gibbs measure and the 
Glauber dynamics and then discuss some basic results on convergence to the plus phase, 
together with properties of the sets il^ introduced above. 

2.1. Finite Gibbs measures on the 6-ary tree. We denote by d{x,y) the tree distance 
between two vertices x, y G T''. If r is the root of the tree, we write d{x) = d{x, r) for 
the depth of x. When A is a subset of vertices of we set d{x^ A) = inf^gA d{x, z). The 
boundary of A, dA, is defined as the set of vertices x such that d{x, A) = \. E{A) denotes 
the set of T^-edges (x, y) with x,y & A. 

The Ising spin configurations space is the set Q, = {—1, +1}^* and its elements will be 
denoted by Greek letters a, r/, ^ etc. The set VL is equipped with the standard cr-algebra 
F generated by the variables {(yx}x&Tb- ^'^^ ^"^Y finite subset ^ C T'' and any rj ^ Vl, we 
denote by /x^ the Gibbs distribution over conditioned on the configuration outside A 
being ry: i.e., if cr G agrees with r] outside A then 



where (i is the inverse temperature and h the external field. We define fi\{a-) = oth- 
erwise. If the boundary configuration r] is identically equal to +1 (—1) we will denote 
the corresponding conditional Gibbs distribution by //^ (/^^)- Whenever the set A will 
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coincide with the finite subtree = {x G : d{x) ^ i} we will abbreviate the symbol 
Ti in the Gibbs measure with £, i.e. fj,^ stands for /ij.^. 

For a bounded measurable function / : 17 ^ M we denote by fj\{f) = Z^^en /^a('^)/('^) 
the expectation of / w.r.t. the distribution Analogously, for any X € n\{X) := 
y^i'^x) where Ix is the characteristic function of the event X. We will write Var^»7^(/) 
or \a.v\{f) for the variance ) - fJ^if? and (for / ^ 0) Ent^^(/) or Ent^(/) for the 
entropy /u^(/log/) - log^^(/) w.rt. Note that Var^(/) = iff, conditioned on 

the configuration outside A being t], f does not depend on the configuration inside A. The 
same holds for Ent^(/). We shall use the symbol for the map r/ Similarly, Var^ 

and Ent^ stand for rj \ai\ and rj Ent^. 

A probability measure /i on will be called a Gibbs measure for the Ising model 

with parameters h) if 

^l{^XA{X)) = fi{X), for all X G and all finite sets AcT'' . 

In this work a crucial role will be played by the following monotonicity property of the 
Gibbs measures (and of the Glauber dynamics, see below) known as attractivity. For any 
increasing bounded measurable function /: 

(i) for any A cT^ the map r/ fJ^if) is increasing; (2.1) 

{ii) i4}{f) ^ ^\{f) whenever AdB. (2.2) 

Recall that the "plus phase" ^+ is obtained as the weak limit as £ ^ oo of /i^. Existence of 
this limit follows from the monotonicity properties described above. Similarly one defines 
the "minus phase" It turns out that any (infinite volume) Gibbs measure ^ satisfies 

/i^ ^ ^ ^ 

2.2. The Heat Bath dynamics on finite trees. For any finite subset ^4 C T'' and any r G 
we define the Heat Bath Glauber dynamics in A with boundary condition (b.c.) r (see e.g 
||18'|) as the continuous time Markov chain on il^ := {a G { — 1, 1}^^^^ : = r on dA} 
with generator 

(£y)(a) = j;c,(^7)[/(a-')-/M, 'y^^A. (2.3) 
where {(T^)y = ay for ally ^ x and {(J^)x = —^x and 



""y ^y. d{x,y)=i yj 



1 + Wx{a) ' 

In analogy with the infinite volume case discussed in the introduction the chain started 
from ^ will be denoted by {erf '"^'^jt ^ o- If ^ = we will simply write erf 

It is well known that there is a global pathwise coupling among the processes 

{((TP'^)t>o, A C T^ G n} such that, for any A c 5 C T^ any C and any 
r < r': 

af'^'- ^ af'^'^ ^ af'^'+ Vt ^ (2.4) 

It is a well-known (and easily checked) fact that, for any finite A c and any r, the 
Glauber dynamics in A with b.c. r is ergodic and reversible w.r.t. the Gibbs distribution /i^, 
i.e. for any function / 

lim e*^A/ = MA(/)- 

t—foa 
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The rate at which the above convergence takes place is often measured using two concepts 
from functional analysis: the spectral gap and the logarithmic Sobolev constant. We now 
describe these two quantities for a generic (finite or infinite volume) Gibbs measure fi. 

For a local function f : U ^ R define the Dirichlet form of / associated to the Glauber 
dynamics with reversible measure /i by 

^A<(/) = ^E^(^4/(^')-/(^)]') =E/^(Var|,.|(/)). (2.5) 

X X 

(The l.h.s. here is the general definition for any choice of the flip rates Cx', the last equality 
holds when specializing to the case of the heat-bath dynamics.) The spectral gap Cgap(/u) 
and the logarithmic Sobolev constant CsohifJ-) of fi are then defined by 



■-gap 



/ Var^(/) f^o Ent^(/) 



where the infimum in each case is over non-constant functions /. 

The spectral gap Cgap(/i) measures the rate of exponential decay as t ^ oo of the vari- 
ance w.r.t. /i, i.e. Cgap(/i) is the (largest) constant such that for any / 

^^T^{Ptf) ^ e-2*^^-p(^) Var^(/) , (2.7) 

where Pt denotes the semigroup associated to the Dirichlet form V^{f). The log-Sobolev 
constant Csob(Ai) is related to the following hypercontractivity estimate (see e.g. II24II ): 
setting qt -.= 1 + e^^=°b(/i)t ^ve have, for any function / and any t ^ 

ll^'JII,.,;.^ ||/||2,M, (2.8) 

where stands for the L^'-norm of / w.rt. ^. 

If /i is s. finite volume Gibbs measure (i.e. ^ = /i^) then both Cgap(^) and Csob(A') ^re 
always strictly positive (possibly depending on A, r). The striking result of [19 ] is that the 
same is true for any choice of the parameters (/?, h) if ^ = ^jl^ is the infinite volume plus 
phase. More precisely one has 

inf Cgap(/u|) > 0, inf Csob{^4) > ■ (2-9) 

Such a result does not imply however any ergodicity statement for the infinite volume 
Glauber dynamics. Simple monotonicity considerations show in fact that for any increas- 
ing local function / and any t ^ 0: 

Ptf{-)^^^-{f)^^^^{f)^Ptf{+), 

i.e. non-ergodicity whenever /i^ / /^+, that is /3 > /3o and \h\ < hc{(3). 

2.3. Convergence to the plus phase: preliminary results. A first important step in the 
proof of Theorem ll.2l is to show that convergence to /i+ occurs when we start from all + 
spins. Recall that pt.xiv) = ^((^tx) stands for the expectation at time t under the infinite- 
volume dynamics started in ij. 

Lemma 2.1. For all 6, (3, h there exist 6 > such that the following holds. For all x G 
there exists to{x) < oo such that ift^ to{x) then 

^ PtA+) - /^^K) ^ exp(-(yt) (2.10) 
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Proof. The left inequality is a direct consequence of monotonicity (see (I2.4D ') and the fact 
that fi^{pt,x) = /^"^(c^x) for all x since /i'^ is an invariant measure. We now prove the right 
inequality. For simplicity we only analyze the case of the root x = r (the general case 
requires no modifications in the argument.) Fix a length scale £ and observe that for any rj 
monotonicity implies pt.riv) ^ Pt'^iv)' with the latter denoting expectation of crj'^''*" (the 
spin at the root at time t for the dynamics in Ti with + b.c. at the leaves of Ti and initial 
condition r/). Let also cpf^ denote the function t] Pt'^iv) — pti^'r) so that 

Setting qt:=l + e^^sobC/^l)*^ the estimate (12. 8D yields 

Let ^/(+) denote the probability of having all + spins in Ti under p^. Then there exists 

Ci < oo such that ^ e^'~^^^\ Moreover by fl2.9D we know that qt > e^^* for some 

positive C2 independent of £. Then 

^ exp{C,b'e-''')Ulfj,^^+. (2.11) 
Set now £ = c^t with cs > small enough. Using (12. 7D and fl2.9D we therefore arrive at 

(pi'^i+) ^ exp (Ci6^e-"2t)g-cgap(/.+)t ^ ^-c^t 

for a suitable constant C4 > and t sufficiently large. Now the claim (12.101) follows from 
the fact (see e.g. [19]) that in the + phase the influence of plus boundary conditions 
decays exponentially fast at any temperature: there exists C5 > such that 

□ 

The previous result allows to show that the set is increasing, i.e. its indicator function 
is increasing. 

Corollary 2.2. For any a € (0, 1) the event is increasing. 

Proof. We need to show that for any pair {ri',ri) with r/' ^ rj and r] G fla, also the first 
component r/' belongs to $7^. To prove the claim we observe that, for any x G and any 
t > 0, monotonicity implies 

pt,x{'n) - p^{(Jx) ^ ptAv') - p^ic^x) ^ pt,x{+) - /^^(-^x) (2.12) 

The l.h.s. of (12.121) is ^ — e^*" for any large enough time t because G il^ by assumption. 
The r.h.s. is instead bounded via Lemma IZTl above . □ 

Another consequence of Lemma IZD is the following 

Corollary 2.3. For any rj ^VL^ the law of the process converges weakly to p^ as t ^ 00. 

Proof Observe first that by the global coupling, for any x G and r/ G 17 we have 
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Let / be a function on ft depending only on the spins in a finite set A cT^ and let r/ e Oq. 
Then, using the invariance of i.e. fi^Pt{f) = J dfi'^{^,){Ptf){£,) = for all t large 

enough depending on A, we have 



lEl/Ca?)) - ^+(/)| ^ |E(/(a?) - /(a+))| + | J d^+(OlE(/K+) - /(4))l 

^ 2 ll/IU ^ [ Fia^x) / atix)) + / / ^*+(^)) " 

= ll/IU E [ 2a.-(+) - ax(r?) - t^-'iax) ] ^ \\f\\oo\A\ [ 2e-^* + e-*" ] . 

Therefore K{f{a^)) fJ-~^{f) for every bounded local function and the weak convergence 
6jjPt iJi^ follows. □ 

Finally, the following generalization ofLemma l2.1l will also be useful. Let us define the 
set ill 5, for (5 > 0, as the set of 77 G such that (11.31) above holds with the stretched 
exponential exp(— t") replaced by the true exponential exp(— Lemma IZTl then says 
that + G $7i 5 for some J > 0. 

Lemma 2.4. For every b, P, h, there exists 6 > such that u{ni^s) = ^for any v ^ /i"*". 

Vroof. Since v ^ /i+ we have iy{pt,x) ^ l-i-'^{(^x) for all t ^ 0. From Lemma IzTl we then 
infer 

u{\pt,x - ^ e--^*/^) < e^*/2 ^(|^^^^ _ ^+(^^)|2) 

Therefore, the Borel-Cantelli lemma implies that there exists e = e{6,b) > and a subset 
ilo C of z^-fuU measure such that for all r/ G JIq) all integers j large enough and all 
X G Tgj (the tree of depth £ = ej), 

|Pj,x(^?)-/U+(a,)Ke-'^^V4 (2.13) 

To prove the lemma we will establish a bound of the type (I2.13D on \pt,r{'n) ~ M^(<7r)|, i-e. 
at the root x = r, but for all times t large enough and not just integer ones. The case of 
general x is obtained by straightforward modifications. We simply write pt for pt^r- Then, 
if [t] is the integer part of t: 

Ptiv) - ^^'^i'^r) = P[t\ir]) - fJ-'^icTr) + dsPsgir]), g := Car (2.14) 

Jlt\ 

For s ^ [tj the Markov property yields Psg{ri) = PitiP.-ingiv) = E([P,_Ltj5](f^Ltj ))• On 
the other hand standard arguments (the so-called "finite speed of propagation" estimate) 
based on tail estimates for the mean one Poisson process (see e.g. [18]) show that 

sup \Pug{0 - PugiOl ^CiV] le. - Q e-""^'^"^ (2.15) 

<u<l 
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for constants Ci , C2 with the property that we can take C2 as large as we wish provided 
Ci is large accordingly. Therefore by the global coupling and the invariance of 

\Psg{r,)-lx+{g)\ = \ j dM+(e)E([P,_Ltj5](^^f,j ) - [^-L^j^K^J ))l 

To handle the last term we add and subtract ctj^j ^ so that by monotonicity 

/ - 4j,J) ^ 2F(.+ = 1) - n^l^^^ = 1) - / rf^+(0F(4j_, = 1) 

Fix j = [t] . When x G T^j we use (I2.13D for the first term above. If e is sufficiently small 
the argument of Lemma ITT) also yields 

P[t\,xi+)- f'^i'^cc) ^ e-^'\ 
for some 5i > 0, uniformly in x G T^^i^jj . In conclusion, for a suitable ^2 > we have 

\P,g{ri)-fi+{g)\^e-'^' (2.16) 

for all sufficiently large t. The desired estimate now follows from (I2.16D and (I2.14D . since 
f^~^{g) = by invariance of □ 

3. Proof of Theorem Ii.2I 

We will provide a unified proof of the three statements in Theorem ll.2l In order to be 
able to do so we need some preliminary observations. The first is that by the monotonicity 
of the events rCorollarv l2.2D . statement b) in Theorem ll.2l is equivalent to 

b *) For every p > \ there exist bo and Po such that for b ^ bo, (3 ^ Po and h = 0, we have 

Fp{Qa) = I for some a = a{ j3,b) > 0. 

Similarly, performing a global spin-flip, statement c) in Theorem ll.2l can be rephrased as 

c*) For every p > there exist bo and Pq such that for b ^ bo, P ^ Po and h = +hc{P, b), we 
have ¥p{VLct) = Ifor some a = a{P, b) > 0. 

The last observation is that we may replace statement a) in Theorem ll.2l with 

a*) For every a > 0, b ^ 2, there exist p < 1 and Po > such that for all /3 ^ /3o and 

h ^ — hc{P, b) + awe have ¥p{Qa) = h for some a = a{P, h, b) > 0. 

In other words, we are taking P large enough. To see why this is not restrictive recall that 
by an obvious domination argument one has Pp ^ ^+ if 

g(6+l+h)/3 I p 



e(fe+m)/3 + e-(6+i+/i)/3' •• ^--2{b + l + h) ^^l-p'' 

Lemma |Z4I therefore implies that fp{Q.a) = 1 for all a < 1 if p ^ Pi3,h- We then achieve 
the result of Theorem ^2] a) from a"") above by a suitable tuning of the parameter p. 
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Figure 2. Free vertices (o) and obstacles (•) in a given realization uj on 
the binary tree. 

3.1. Main argument. As the convergence result of Lemma l2.1l makes clear, to prove The- 
orem ^2] we need to lower bound the quantity pt,x{'n) ~ in the three statements 
a*),b*) and c*) emphasized above. We shall focus only on the case x = r, since the case 
of arbitrary x is obtained with essentially no modification. Setting ptit]) := pt,r{i]) what 
we want is a bound of the form 

pt{T]) - p+{ar) > - e^*'' , (3.1) 

for all t ^ to{ri), Pp-almost all 77. As far as this section goes, we shall not distinguish the 
specific setting (a*,b* or c*), since all we do here works for the three cases without any 
difference. What does depend on the setting are some key estimates that will be proved 
in the next two sections. The latter have been emphasized as separate claims in the text 
(see Claims 1 to 4 below) . 

In order to describe the main idea behind the lower bound (13. IB . we must first in- 
troduce the notion of the Ising model and the associated Glauber dynamics in a random 
environment of obstacles. Realizations of the environment are described by elements uj of 
n. We say that a vertex x G T'' is an obstacle if uJx = —1, and that x is free if lox = +1- We 
call T{uj) the largest connected component of the set of free vertices containing the root. 
Note that T{uj) = if the root r is itself an obstacle. By construction, all vertices in dT(uj) 
are obstacles. We will be mostly concerned with the case where uj is picked according to 
the product Bernoulli measure Pp, i.e. when each vertex is free with probability p, inde- 
pendently of all others. In this case, Pp(T(w) is infinite) is positive as soon as p > 1/6 and 
tends to 1 as p X 1 for fixed 6, or as 6 ^00 for fixed p, see e.g. [23]. 

Given a realization of obstacles uj, the Ising model among obstacles is defined as before 
by replacing the tree with the random tree T{uj) and the configuration space $7 with 
the space 

:= {r G : = -1 , Vx ^ T{uj)} . 
Given a finite subset A <zT^ and r G we denote by the Gibbs measure p\p^rp{uj)- 
We also write /ij^ for the Gibbs measure Pxi{uj)' where we use the notation T(^{tjj) := 
n T{ijj). From this definition we see that obstacles act as a "minus" boundary condition. 
The maximal allowed configuration t ^ is such that r = +1 in T{uj) and, when no 
confusion arises, it will be always denoted by We will write for the Gibbs measure 
obtained as weak limit of p^^ as ^ ^ 00 (this is a finite volume Gibbs measure with — b.c. 
if r(w) is finite). 

Similar notations apply to the Glauber dynamics. Given a realization of obstacles uj, 
^ C and r G B^, we will write af^''^ for the Glauber dynamics in ^ n T{uj) with 
boundary condition r started from the restriction to T{uj) of the configuration ^ G SI. If 
yl = we simply write crf^'^. When A = T{uj) the boundary condition is necessarily "— " 
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and we will only write a^^. In the same way we will use pt, w{0 the expected value 
of a^^ at the root and pfljiO for the expected value of cr^'^'^ at the root. Monotonicity 
implies that, for any ^ , cr^ ^ < . In particular, 

Pt,U+) ^ Pti^) ^ i^O- (3.2) 

We now turn to our main argument. Fix a length scale £, to be related later on to the time 
t, a configuration 77 G and define the associated realization uj = uj{rj,i) of obstacles by 
the rule: 

= "<^><^ (3.3, 

I r]x otherwise 

Clearly crj'^ ^ crj' so that 

Ptiv) - P'^{(^r) > [Pt,cj{v) - P^i^^r)] - [p'^{(Tr) " pi{(^r)] ■ (3.4) 

If now L = P , J > 1, is another length scale, monotonicity shows that if we impose + b.c. 
on the leaves of Tl{oj) we may estimate 

> [Pt'^iri) - pUM] - [pt'J'iv) - Pt,M)] - [P^i<^r) - Piiar)] . (3.5) 

Notice that in the above formula the role of the Bernoulli configuration rj is twofold: it 
enters as the starting configuration in the first two terms but it also defines the random 
realization of obstacles uj. 





+ + + + + + + 

Figure 3. Random obstacles below level £: Infinite tree (left) and finite 
tree with + boundary condition below level L = P, j > 1 (right). 

Most of the statements that will be proved below on the rh.s of (13. 5D concern properties 
which hold almost surely with respect to the starting configuration ij (and therefore wrt. 
uj) picked according to the Bernoulli measure Pp. To simplify the exposition, we shall 
adopt the following convention: given some statements ££, £ G N, we say that £e holds 
Pp-a.s. for £ sufficiently large whenever r] e £e for all £ ^ £o{'r]), for Pp-a.a. G $7, or 
in other words, fp{£i eventually) = 1. We are now in a position to explain how we will 
bound the three terms in the rh.s of (13. 5D . 

Estimate on [p'^ (ar) — p^ (err)]- Bounding the third term in (13.51) is a purely static problem 
which on the tree can be solved via a suitable recursion. In section[S]we prove 

Claim 1. 

(3.6) 

Pp-a.s. for £ sufficiently large. 
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Estimate on [ptlj~{ri) - /i^^((Tr)]. The first term in (I3.5D is related to the speed of relax- 
ation to equilibrium in the finite tree Tl{uj) with plus b.c. Here we need the following 
bound on the logarithmic Sobolev constant Csob(^L 

Claim 2. There exists ( < oo independent of £ such that 

holds ¥p-a.s. for any L ^ I sufficiently large. 

We can repeat the computation (12.111) with fij' replaced by fJ-'^^, (pf^ by (f't'J'ir]) := 
Pt^^iv) - Pt,uji^r) and qt:=l + g^^'^^^bC/^t-) to obtain 

\c^f:^^{v)\ ^ exp (Ci6^e-^-('^t J*) ll^^+JI^,^^^ ■ (3.7) 

for some constant Ci < oo. Assuming Claim 2 above, using Cgap ^ 2csob (which is always 
true, see e.g. [24]) we estimate (13. 7D with the help (12. 7D and obtain, for L = P, 

Ipt^iv) - Ptj^r)\ = < exp (Ci6^e~*/^^^)e-*/^^^ (3.8) 

Pp-a.s. for £ sufficiently large. 

Estimate on [pf ^ {rj) — pt^ui (??)]• The control of the second term in (13.51) is a true dynamical 
question and it involves proving that the two processes u^^^ and cr'f '^'^ remain identical 
at the root up to time t with large probability. This is achieved via a coupling argument 
together with some equilibrium estimates. The final bound will be of the form 

Pt^{^)-PtM^te-''' (3.9) 

Pp-a.s. for I sufficiently large, provided that 7 > C + 1, where C, is the constant appearing 
in (13. 8D . The argument goes as follows. 

To keep the notation to a minimum, we will abbreviate the two processes and cTt'^'^ 
with S^j and respectively. Set 

A:={x€ T{uj) : d{x) = 21}, A := {x G T{uj) : h ^ d{x) ^ ^i} . 

By the global coupling 



^ 2] 



3s^t,3xeA : CL/e 



2 

s.x 



(3.10) 



Define Aj, j = 1, 2, . . . , [tj, as the event 

Aj = {3xeA : 
The r.h.s. in (I3.10D is then estimated from above by 

Ltj L*J+i 

^F(^j-) + (3.11) 

where 

Bj := {3s G [j - l,i] , 3x G A : + ^s,x } • 
The probability of the event D Bj is estimated by a standard argument: the event 
^j-i ^ implies that a discrepancy between and located outside A, reaches in 
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a time smaller than 1 a point x G A. Since there are at most possible (self-avoiding) 
paths from A to (A)^ and since the rates are bounded by one, a simple tail estimate for 
Poisson random variables implies 

F{A'j_^ n Bj) ^ c6t^e-5"°gW^) ^ e'^^ (3.12) 

for a suitable constant c and all sufficiently large £. Similarly, if we look at the event Aj, 
we are requiring that at least one of the discrepancies at time in level L travels up to 
level |£ in a time less than j. Therefore 

for some c < oo and for all e = e{b, c) sufficiently small. We are therefore left with the 
estimate of P(Aj) for j ^ eL. The argument for this case goes as follows. 

Fix a point x G A and recall that x is at some level between |£ and |£. Let = r{x, £) 
be the ancestor of x at level £ and let Tr^{uj) be the subtree of T{uj) rooted at r^. and 
containing all descendants of r^. Let also T^^ /i e N, be the finite subtree of Tr^{u}) 
obtained by considering only the first h levels of T^^ (cj) (so that Tj.^^h{uj) consists of rx and 
all its descendants in T{ijj) lying between level £ and £ + h). When h = 21, Tr^^2(.{^) is the 
tree between levels £ and M, so that A c Tr^,2^(^) and c?(A, {Tr^,2t{'^)Y) ^ f ■ Call t'.^'J^ 




3^ 



Figure 4. The vertex x e A and the associated tree Tr^^2e- 

(f^'J^) the Ising-Gibbs measure on Tr^^hi^^) with + (— ) b.c. above the root and + b.c. 
below the leaves at level £ + h and let i^rjj^ = lim/i_^oo ^^'t- 

Finally, we denote by the Glauber dynamics evolving in Tr^^2e{^) with + bound- 
ary conditions both above the root and below the leaves of T,.^ 2^('^) and with initial 
configuration -q. Notice that in fact starts from all pluses because, by construction, 
= +lVy er(a;), ^ 

With the above notation and using monotonicity we can now write 

^ ¥[Clx = +1] - (a, = +1) . (3.14) 
The r.h.s. in (I3.14D is then decomposed into the sum of two terms: 

^[^Ix = +1] - = +1) > (3-15) 

and 

= +1) - ^rjU^. = +1) • (3.16) 

In order to bound the term in (13.151) we would like to argue as in (13. 7t and therefore we 
need the following: 
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Claim 3. There exists ( < oo such that 

min Csoh{i^r'te) ^ ^ 
xeA 

holds ¥p-a.s. for I sufficiently large. 

The argument used in fl3.7D now yields that the first term (13.151) satisfies 

= +1] - = +1) < exp {cb^'e-^/'') exp ( - . (3.17) 

Therefore, if 7 > C + 1> using j ^ eP we may write 

n^l. = +1) - = +1) ^ b-^'e-^' (3.18) 

for £ large enough (independent of x). In conclusion, for any 7 > C + 1, 

E {ml- = +1] - = +1)} ^ (3.19) 

xeA 

Pp-a.s. for ^ sufficiently large. 

As far as the term (13.161) is concerned we will establish: 

Claim 4. 

E = +1) - ^rjU'^x = +1)} < e-^*^ , 

xeA 

Pp-a.s. /or £ sufficiently large. 

Collecting ( I3.19D and Claim 4, we have thus shown that P(Aj) < 2e"^^, j ^ eL = ef . 
Together with (13.121) and (13.131) this completes the proof of (|^?9t . 

Conclusion. In conclusion, from (13. 51) . using (13.61) . (13.81) and fl3.9D we have the bound 

p,{r^) - ^ - + exp (C7i6^e-*/^^')e-*/^^' + t e'^^) . 

If e.g. ^ = t^ with Q > such that a7(l + C) < 1 and 7 > ( + 1, then 

pt{r])- fi+{ar) ^ -e-*" , (3.20) 

Pp-a.s. for t sufficiently large. Therefore the three steps above are sufficient to end the 
proof of Theorem 11.21 Note that the coefficient a depends on the various parameters 
(6, /3, h etc.) only via the constant C, coming from the logarithmic Sobolev inequality in 
Claim 2. 

4. Recursive analysis among obstacles 

In this section we establish a number of key estimates for the Ising Gibbs measure 
among obstacles. Once these results are established it will be rather easy to prove Claims 
1 to 4 (see next section) . 

As in the previous section l<j G $7 will denote a random realization of the obstacle- 
environment and ^+ the associated Ising plus phase. We emphasize however that here, 
contrary to (13. 3D . w is picked according to the product Bernoulli measure Pp on the whole 
tree T^, i.e. each vertex x G T'' is free with probability p independently of all others. 



16 



E CAPUTO AND F. MARTINELLI 



4.1. Coupling coefficients and path weights. As in the homogeneous case treated in 
II19IL the analysis of equilibrium properties is reduced to the study of certain coupling 
coefficients. For a given oj we define the ratio 

We agree that R{uj) = oo ii uj^ = —1. For every z G T'' we set Rz{oj) := R{0zOj), where 
9z denotes the shift induced by the natural group action on the tree: {6zOj)x = ojz+x- If 
ujz = {6zUj)r = +1 and zi, . . . , z;, denote the children of z G T{oj), one has the following 
easily checked recursive relation (see e.g. [19]): 

b 

Rz{u:)=e^J{Fp{Rz,{u:)), (4.2) 

k=l 

where, from now on, we use the following notation 

F,{a) = ^, e:=e-'^ (4.3) 
1 + ea 

To illustrate the use of the variable R defined in fl4.1D . consider a vertex z G T{uj) together 
with one of its ancestors y and denote by /^^'^ (resp. /i^'~) the measure /i+ conditioned to 
have ay = +1 (resp. ay = —1). Suppose we want to compute the total variation distance 
between the marginals at the vertex z, which we denote by — fiHj" \\z- Since the spin 
at z can take only two values, the latter equals ^M^^{az = +1) — i-i^]~{<yz = +1). If y is the 
parent of z, using ^lj^{az = +1) = {e^Rz{uj) + we see that 

\\^^'J^-^i'i]~\\z = Kp{Rz{u)), (4.4) 

where the function Kp : [0, oo) [0, 1] is defined by 

Kp{a) = - -z^—rr . (4.5) 

For every £ G N we define the set of descendants of y at depth t. 

Dy/{io) = {x G T{u>) descendant of y : d{y,x) = £} (4.6) 

To compute the total variation distance — H^: for some x G Dy^e{uj) we may pro- 
ceed as follows. Let zi, . . . , Z£ = x be the vertices along the path from y to x. We couple 
the measures ^xt]^ ,^xt]~ recursively in such a way that, for every i < i, given that the 
corresponding configurations coincide at z^ then they coincide at Zj+i with probability 1, 
while given that there is disagreement at Zj then disagreement persists at Zj+i with prob- 
ability — /^S' lUi+i = Ki3{Rz-^j^{uj)). In this way the probability of a disagreement 
percolating down the tree from y to x equals 

e 

M'+-fil'-h = llKf,{RzM) (4.7) 

1=1 

Moreover, if |cj — a'\y/ denotes the Hamming distance (counting the number of disagree- 
ments) between a and a' restricted to the set Dy^i{uj), the above argument implies that 
we can find a coupling v^^ of ^^i'^, ixl]^ such that the expected value of \a — a'\y^i satisfies 

lyUW - c7'\y,e) ^ WU^y,.), (4.8) 
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where we introduced the path Ty^^ between y and x, consisting of the sites zi, . . . , ze = x 
as above, and the associated weight 

e 

W{Ty^,,u;) = J]K^(i?,,(c^)) . (4.9) 

1=1 

The rest of this section is concerned with estimates showing that, in a suitable sense, R{uj) 
and Wuj are small with large probability. 

4.2. Estimates on R. We write Fp for the probability Fp conditioned to have tOr = +1. 
We want an estimate of the type 

Fp{R^e) i^6, (4.10) 

where e = e^^'^ and 5 is a small parameter We start with the setting of statement a* in 
the proof of Theorem ll.2l 

Lemma 4.1. For any 5 > 0, a > 0, b 2, there exist po < 1 and /Sq < oo such that 14. J OP 
holds for all p ^ Po, P ^ Po arid — hdP) + a. 

Proof. For any integer i we define 

«V) = ^fc-^. (4.11) 

f^U'^r = +1) 

Since fxj'^ we have R^ ^ R, £ ^ oo, Pp-a.s. Moreover, monotonicity implies 

R^{(jj) ^ R^^'^{uj), so that the convergence is monotone. Then it is sufficient to establish 
fl4.1QD for R^ in place of R, uniformly in L We will give the proof only in the case 6 = 2, 
since all the estimates below are easily adapted to the case of larger values of h. Recall 
that in general (see e.g. [9]) one has hc{P) = (6 — 1) + 0(/3^^), so that, replacing a with 
2a and taking /? sufficiently large, we can assume /i ^ — 1 + a without loss of generality. 
Let us define the probabilities 

qf^ = Fp (r^ > 2-2*^ei-*^") , A; = 0, 1, 2, . . . (4.12) 



Let now zi, Z2 denote the two children of the root and observe that the corresponding 
ratios Ri, i = 1,2 are i.i.d. random variables with the same distribution as R^~^. On the 
event {cor = +1} the basic relation (I4.2t applies and we have 

R^ e^-^ Fp{Ri) FpiR2) . (4.13) 

Using the uniform bound Fjj ^ e~^, we see that, in particular, R^ ^ e"^^e^^. Therefore 

q^'^^ = for every £ ^ 1 as soon as A; > fco := L| ~ IJ ^^id /? is large enough. Suppose now 
Ri ^ e. Then by (l4!T3l) . using Ffs^Ri) ^ 2e we have 

R^ 2e''Fp{R2 



12 J 

-1- 



Since Ffs{t) ^ e+t, the event R e forces R2 ^ for large enough /?. Considering also 
the event {w^-^ = —l}VJ{ujz2 = ~1} and the event {wz^^ = +l,i?i > £}r\{ujz2 = +1,-^2 > e} 
we may then use (14.121) to write 

gf)^2(l-p) + (,(V + 2,5i\. 

The same reasoning as above actually shows that for any k one has 

,f ^2(l-p) + (,f_\)2 + 2£t^). (4.14) 
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From the monotonicity in I of we see that qf}_^ ^ q^^ for any k and I. Therefore a 
simple iteration of (I4.14D gives that 



i-i 

,(i) 

m=0 



^ ^2'"{2(l-p) + (gf_\)2} + 2^ 



for any j = 1, 2 . . . . When j = fco + 1, g^i^ = and we have the recursive estimate 

^2^0+1 {2(1- + (4.15) 

This implies that for every 5 > we can choose po < 1 and /3o < oo such that qy' ^ 5, 
for every ^ ^ 1, p ^ po and /3 ^ /3o- To see this simply observe that when i = \ the "+" 
boundary condition imposes i?j = on every child Zi such that a;^. = +1 and therefore 

qi' ^ 2(1 — p), which can be made arbitrarily small. Thus, climbing up the tree with the 

relation (I4.15D . we see that sup^ > iQi ^ (5 as soon as e.g. p ^ 1 — 1/ (2^^^^^^"). D 

We turn to the setting of statement b " in the proof of Theorem ll.21 

Lemma 4.2. For any 6 > 0, p > ^, there exist &o S N, /3o < oo and c > such that H4. 1 OH 
holds for all p ^ Pq, b ^ bo and h = 0, with 6 = e'^K 

Proof. Recall the definition fl4.11D of R^{uj). Let zi,. . . ,Zb denote the children of the root 
r and let m{uj) stand for the number of obstacles among them: m(w) = ^2^1=1 =-i}- 
Since p > ^, from standard large deviation estimates for the binomial distribution there 
exist positive numbers ai , 02 > and 60 £ N such that such that 

m ^ ( ^ - ai ) b] ^ e-"2'' ^ (4.16) 



2 

for all b ^ bf). Suppose now that uJr = +1 and m < (^ — ai)b, i.e. the root has at least 
I + ai6 free children. Suppose only one of these free children, say z, is such that the 
associated ratio Rz satisfies Rz ^ e. In this case (14.21) yields 

R ^ e~'^i2ef-'^-^F{Rz) ^ 2^e2aib-i^^(^^) _ (4_;^7) 

Since Fp ^ e^^ it is clear that we can take 60, /3o so large that in the above situation it 
is impossible to have R ^ e for all 6 ^ 60 and P ^ Pq. The above discussion says, in 
particular, that if m < — ai)b and R ^ e, then there must be at least 2 children of r 

with ratio smaller than e. Thus, recalling the definition of the probabilities qf^ (I4.12t . we 
obtain 

+ ^ (')(,(«.)". (4.18) 

n=2 ^ ^ 

Because of the + boundary condition at level £, the argument of (I4.17t gives qf'^ ^ e""^^. 
The claim then follows by induction: Suppose q^f]-^ ^ 6 with 6 := e""^''/^. Then (I4.18D 
implies ^ ^ e-'^^fe + (1 + ,5)^ - (1 + 5b) ^ e-'^^fc + 1 52^2^ therefore ^ ^ 5 for 6 
suitably large. □ 

Finally, for the statement c* in the proof of Theorem ll.2l we need the following 

Lemma 4.3. For any 5 > 0, p > 0, there exist bo £ N, Po < 00 and c > such that 14. JOD 
holds for all P^ Po,b^bo and h = +hc{P), with 5 = e'"*. 
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Proof. As in the previous proof we denote by m{Lu) the number of obstacles among the 
children of the root. Since p > 0, there exist positive numbers oi, 02 > and bo such 
that 



for all b ^ 6o- We recall that hc{f3) = (6 — 1) + 0(/3~^). In particular, we may assume 
without loss of generality that the magnetic field satisfies h ^ 6 — 2. If zi, . . . , denote 
the children of the root, by (14. 2D we then have 



Reasoning as in (I4.17D we see that if = +1 and m{uj) < (1 — ai)6, then we must have 
more than one free children with ratio smaller than e in order to produce the event e. 
It follows that we may estimate the probabilities qf^ exactly as in (I4.18D . When i = 1 the 
+ boundary condition implies R ^ e^^'^e^"^e^~"^. Therefore on the event m{uj) < (1 — ai)b 
it is impossible (for suitably large b,(3) to have R ^ e. This gives q^'^ ^ e^^^f*. As in the 
proof of Lemma lT2l the desired result now follows by induction. □ 

4.3. Estimates on W. We turn to an estimate on the weight W introduced in (I4.9D . Re- 
call the definition (I4.6D of the set Dy i{uj). Below we simply write = D^{uj) when y 
coincides with the root r. We also write W{x) := W{Tr^x^^)i, for any x G T^. We look for 
an estimate of the form: There exists to > such that 



for every t ^t^, \. The value of to will depend on the parameters a, b, 13 in case a*, on 
6, 13 and p in cases h*,c*. 

We start with the setting of statement a*. 

Lemma 4.4. For any a > 0, b ^ 2, there exist /3o < 00 and po < 1, such that H4.20\l holds 

for any (3 Po, P > Po, h ^ - hdP) + a . 

Proof The main difficulty in proving (I4.20D is the non-independence of the random vari- 
ables {Kf^{Rz{uj))}z^ry X entering in the definition fl4.9D of W^j. However, thanks to the 
tree structure of our graph, it is possible to introduce a modified weight Waj{Ty^x) ■= 
TlzeTy X in such a way that: 

• the random variables {ipz{uj)}z&ry ^ are independent; 

• for each (J, WiJ^y^^.u)) < W{Yy^^,u)). 

We now describe how we construct the modified weights. 

To begin with, we fix some notation: dBi stands for the (deterministic) set of vertices x 
such that d{x) = £. We also use for the unique path from the root to x. To simplify the 
notation we define Wuj{x) = if x G dB^ \ Di{uj). Next, for every x G dBi and for every 
vertex z G F^^ we denote by the set of all children y of z such that y ^ T^- Clearly, 
I A^l =6 — 1. We say that z G T^; is regular if Ry ^ e for every y G A^. 

Let now u > be a small parameter to be fixed later and suppose that zi, Z2, . ■ ■ , z^. are 
consecutive regular sites on Fx, ordered in such a way that d{zj,r) = d(zj_i,r) — 1. Let 
also zo denote the children of zi along r^;. As in the proof of Lemma RTTl we may assume 
h ^ — (6— l) + a without loss of generality. Since zi is regular, using I-F/3I00 ^ 



P(m ^ (1 - ai)6) ^ e"''^ 



b 



(4.19) 



R^e''-^FpiRz,)---Fp{Rz,). 




(4.20) 
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Fj3{x) ^ 2x, for any x ^ e, from (I4.2t we have 
Therefore, by proceeding inductively, for any k we have 



^ ^a-(fe-l)^(fc-l)a-l2{fc-l)6(2e)''-l 



(4.21) 



Suppose now k ^ ko '■= [|J. The above argument shows that R^^ ^ ne, provided 
e ^ eolO)^)?^)- A simple computation gives K{ae) ^ a for every a > 0, so that 

Kp{R,,)^u. 

We shall say that z G F^; is good if z is regular and the number of consecutive regular 
vertices immediately below z along is larger or equal to /cq — 1- Otherwise we say that 
z is bad. 

The estimate (I4.21D therefore implies that Kji{Rz{uj)) ^ u whenever z is good. Since 
Kp ^ 1 and recalling that W(x) = if x is not connected to the root in T(u;) we may write 



U UJz 

1 u:z 
uj. 



+1 , zis good 
+1 , z is bad 
-1 



(4.22) 




(4.23) 



Figure 5. k consecutive regular sites on the path in the case 6 = 4. 

We now claim that there exist Ci < oo such that for every £ ^ 0, for all x G dBf. 

Ep[W{x)] ^ Ci{2uY, 
From (I4.22t we see that 

Wix) ^ n^-"-('^) . (4.24) 
where n^ (w) stands for the number of bad vertices in F;^.. Define now, for every z G F^^, 
Xz{^^) = if z is regular and Xzi^^) = 1 otherwise. Note that, by construction, these are 
i.i.d. Bernoulli random variables. A simple deterministic bound on Ux is given by 



H ^ fco (i+ 5^ X.M) • 



From Lemma im we know that the probability of being irregular, for any given z G F^^, is 
less than 6i := (6 — 1)6 + (6 — 1)(1 —p). Let us choose 6 in Lemma ITTl and p < 1 such that 
6i ^ . We then have 



W{x) 



-ko 



u 



^ U 



e-koi 



(4.25) 
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The claim (I4.23D then follows by taking Ci = u~^°. 

We are ready to prove the exponential moment estimate (I4.20D . For any integer k we 
define 

Mfc = j W{x) j = ^(^i) • ■ ■ ^(^'^) ■ (4.26) 

We claim that 

Ep [Mk] ^ C^k\ , = 1, 2, . . . (4.27) 

for some constant C2 < oo. Note that the result (I4.20t is an immediate consequence of 
(I4.27D since the l.h.s. in (I4.20t is bounded by 

°° 1 
Ep [exp (tMi)] = 1 + ^ ^^P^'^ ^ T3c^ • 
fc=i ■ ^ 

Let Xl, . . . , Xfc e be given as in a generic term in the sum in (I4.26D . These points may 
be ordered by the lexicographic rule to obtain the ordered set xi ^ 5:2 ^ • • • ^ Xk- Call 
Xo and Xk+i the absolute leftmost and, respectively, the absolute rightmost vertex in dB^. 
Below we use [xj_i, Xj) to denote the set of vertices y G dBi such that y is larger or equal 
to Xj-i but strictly less than Xj, with the agreement that, when j = k + 1, the set [xfc, x^+i) 
also includes the end point Xfc+i- With these notations we can write 

k+l 

j=l Xi_,...,x^^dBi y£[xj-i,Xj) 




Figure 6. Schematic picture of the tree Tj. in the case k = A. Here y e [x4, X5). 

Consider now a given y G [xj_i, Xj). Let Tk denote the subtree determined by the union 
of all paths Tr^xi, i = I, ■ ■ ■ ,k. Let d{y, Tk) denote the distance from y to and write Zy 
for the closest ancestor of y on Tk (characterized by d{zy,y) = d{y,Tk)). Clearly, we can 
estimate W{y) ^ W{Tzy^y), where W{Tzy,y) = Ilzer^ y Now, by construction, the 

random variable W{Tzy^y) is independent of all the weights W{xi) except for the variables 
ipz where z is either Zy or one of the ko — 1 consecutive vertices just above Zy. Let us call 
Ay this set of vertices. Restricting to the event that xi, . . . , x^ G D(^{tjj) we can estimate 



n m^)r' 



Therefore we have 



k+1 

Ep [Mk+i] ^ u-^^' J2 [^(^1) • • • ^(^k)] Y 12 [^(r-.,: 

Xi_,...,Xk&dBi j = l y£!^Xj-l,Xj) 



(4.29) 
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Clearly, for every integer d there at most vertices y such that d{y, Tk) = d. Therefore, 
choosing u ^ 1/(46), from (I4.23D . for every pair Xj, we have 



ye[xj_i,ij) 



^ 2Ci . (4.30) 



From ( IT291 ) and dTSOb . setting C2 := 2u~''°Ci we obtain 

[A4+i] < C2(A: + 1) Ep [M^] , (4.31) 
so that, for every /c ^ 1 we can estimate Ep [Mk] ^ A;! as claimed in (I4.27D . □ 

We turn to the setting of statement b". 

Lemma 4.5. For any p > ^, there exist /3o < 00 and bo G N, such that H4.20\l holds for any 

P>(3o,b^bo,h = 0. 

Proof. The proof is essentially the same as that of Lemma 14^41 but we have to modify the 
definition of good and bad vertices. Given x e dBe and z ^T^v/e write as before for 
the set of children of z lying outside of the path F^. We write also m{uj) for the number 
of y G such that ujy = —1. As in the proof of Lemma 1431 we may use (I4.16D to 
estimate this quantity: there exist ai, 02 > such that Pp (m ^ (1 — 2ai)6/2) ^ e~°'^'', for 
all sufficiently large b. 

Here the definition of good vertices goes as follows. We say that z is good if m ^ (1 — 
2ai)6/2 and if all the vertices y e such that ujy = +1 satisfy Ry ^ e. Clearly, if z is 
good, from (14. 2t we must have 

In particular, for any u > we find 60 and /3q such that for all b ^ bo and P ^ /3o we have 
Rz ^ ue. As in the proof of Lemma RT?! we therefore have that Kp(Rz) ^ u whenever z 
is good. Now we can define individual weights V'z exactly as in (14 2 2D and, as before, we 
can estimate W{x) ^ W{x). To establish the analog of (I4.23D we simply observe that V'z 
are i.i.d. random variables with the present definition of good vertices. Moreover, from 
Lemma I4T2I we easily infer that 

Pp[2isbad] ^ e"^'' (4.32) 
for some c > 0. Using this we have, see (I4.25D 



Ep 



W{x) ^ {2uY , (4.33) 



as soon as e^^''n^^ ^ 1. The rest of the proof goes now exactly as in Lemma 14^41 The 
estimate (14.291 ) is actually simplified by the fact that we only need to remove the vertex 
Zy, SO that the factor u^'^'^ is now replaced by n~^. In particular, (I4.31D now holds with 
the constant C2 = 2u^^. □ 

It remains to prove (I4.20D in the setting of statement c". 

Lemma 4.6. For any p > 0, there exist /3o < 00 and bo G such that S4.20\l holds for any 

P;? f3o,b^bo, h = +h,{0). 
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Proof. As in the proof of Lemma |43I we call m{uj) the number of obstacles among the 
children in for a given vertex z. We shall use the analog of estimate ( I4.19D . Letting 
ai and a2 be the parameters appearing there, the vertex z is now declared good if m{uj) 
satisfies m{uj) ^ (1 — ai)b and all the vertices y e such that ujy = +1 satisfy Ry ^ e. 

With this definition of good vertices, the bounds of Lemma \43\ now show that 

Pp[2isbad] e'"^ 

for some c > and all sufficiently large b. On the other hand, reasoning as in the proof of 
Lemma 1431 an application of (14 .20 gives that if z is good then 

Given n > we then find feo^/^o such that Kp{Rz) ^ u whenever z is good, as soon as 
b ^ bo and /? ^ /?o. The rest goes exactly as in the proof of Lemma |43] above . □ 

4.4. Poincare and Logarithmic Sobolev inequalities among obstacles. Recall the defi- 
nition (I2.6D of the constants Cgap and Csob- We shall focus here on the case of the measure 
/i = /i^^, i.e. the Gibbs measure with plus boundary condition below a certain level L 
among the obstacle environment uj. We shall write Cgap(i,w) and Csoh{L,uj) for the asso- 
ciated constants. Here uj will be distributed according to Bernoulli(p) measure. 

It is well known that, in general, Cgap ^ 2csob- On the other hand, for trees, a useful 
inequality established in [19] states that Cgap < O(logn)csob) where n is the cardinality 
of the tree. In particular. Theorem 5.7 in [19] in our setting implies that for every b and 
every /3 there exists a constant C < +oo such that for every G $7 and for every L 

Cgap(L,u;) ^ C Lcsoh{L,iL!) (4.34) 

Our main result here is an almost sure polynomial bound on Csob(-^^,'^): There exists a 
constant ( < oo such that 

Csob(^,w) ^L-^ (4.35) 

holds Pp-a.s. for L sufficiently large. Here the constant ( will depend on the parameters 
a, 6, /3 in case a*, on b, /? and p in cases b 

Theorem 4.7. 

a*) For every a > 0, b ^ 2, there exists p < 1 and Po such that 14.351 ) holds for all /? ^ /?o 

and h ^ — hdP, b) + a. 

b") For every p > \, there exist bo £N and Po < oo such that holds for h = 0, b ^ bo, 

c*) For every p > 0, there exist 6o G N and /3o < oo such that 14.351 ) holds for b ^ bo, (3 ^ Po 
and h = +hc{P)- 

Proof We will carry out the proof of the three statements simultaneously. Indeed, the key 
estimate we need is the exponential integrability (I4.20D . which holds in all cases under 
consideration as worked out in Lemma Lemma H3] and Lemma 1431 

Thanks to the deterministic bound (I4.34D it suffices to prove the claim fl4.35D with 
Csoh{L,uj) replaced by Cgs,p{L,uj). We fix a length scale ii much smaller than L. For each 
vertex x £ Tl{uj), let B^/^ C Tl{uj) denote the subtree (or "block") of depth £i - 1 
rooted at x. In this way B^^ii consists of li levels and we understand that if x is A; < 
levels from the bottom of Tl{oj) then B^^i^ has only k levels. In the end we will choose 
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£i = C log L for some sufficiently large constant C. We define the Dirichlet form of the 
so-called "block-dynamics" 

A standard argument relating the spectral gap of the heat-bath dynamics to the spectral 
gap of the block-dynamics (see e.g. [18]) shows that, since there are at most £i blocks 
containing a given vertex x, we have 

Cgap(L, ^ ^ jmin cgap(/i^^^^ )| inf y^l'^^^^^^^ ' C4.36) 

where Cgap(^B ^ ) denotes the spectral gap of the heat bath dynamics on the block -B^/i 
with boundary condition r (and r is assumed to be compatible with the obstacle realiza- 
tion, i.e. T G B^): 



In general trees, according to Theorem 1.4 in [ 1 ], one has a lower bound on Cgap of order 
n^^ uniformly over the boundary condition, where n is the cardinality of the tree and 
C < oo is a constant depending on the parameters b, /3, h. In particular this implies that for 
all w G and for all sufficiently large £i 

- mmCgap(M^^. (4.37) 

Let Cgap(^i, oj) denote the spectral gap of the block-dynamics, i.e. the infimum appear- 
ing in (14.361) . So far we have obtained the deterministic bound 

Cgap (L, a;) ^ 6-2?^^ Cgap(^i,L,a;) • (4.38) 

Next we make a deterministic estimate on Cgap(^i, L, uj). To this end we use the method of 
||19||, combined with the results we obtained in previous subsections. Given r G (0, 1), we 
say that is (£i , r^i )-mixing if for every x G Tl{uj) 



Var^+„ [^^Ui'^x I ) ^ r'^y^r^lJcr,) , (4.39) 

where fi^ ^ {a^ \ (Jd^ ) denotes the conditional expectation of Ux given the values of 
a on Dx/i{ijj), the set of descendants of x at distance li. A simple computation shows 
that fl4.39t is actually equivalent to the variance mixing condition VM(fi, e), with e = r^^, 
introduced in [19]. In particular. Theorem 3.2 in Iil9il implies that 



Cgap(^i,L,w) ^ -, (4.40) 



for ^ Iq, for some finite £q = io{r) as soon as fij^^ is (^i, r^i)-mixing. The conclusion of 
the theorem therefore follows from (I4.38D and (I4.40D if we can prove that //^^ is {£i,r^^)- 
mixing Pp-a.s. for some r < 1, when £i = C log L, with some C < oo, for all sufficiently 

large L. To prove this we observe that, setting ^(cr^) := /ij^ f^L,uji'^x I c^r>,,£,((^)) ka; , we 
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may write 

= 2filJa, = +l)filJa, = -l) [g{+l) - g{-l)] . (4.41) 

Let u denote a coupling of the measures 1 cxa; = +1) and fJ-i^^i-lcTx = —1)- We then 

write 

g{+l)-g{-l) = ^u{T,ri) [^L,^((Tx|Tz5,_,^(a;)) 

with [rr?]^'^ denoting the interpolation between r and t], i.e. the configuration such that, 
using lexicographic order on ^^^^^(u;), {[Tr]]^'^)^ = t^, z < y, and i[T7]]y'^)^ = rj^, z > y, 
while ([Tr7]^'='=)j^ = ±1. Recall now the definition (I4.5D of the function and set 

7 := supK/3(x) = tanh/3 . 

Since [r7/]^'+ and [rr?]^'^ differ only at y G Dx^g-^{uj), reasoning as in fl4.4D and (I4.7D we 
estimate 

I - A^t^K I (N]''-)z>...,h) ^ 27^^ (4.42) 

From (I4.8D we then obtain 

^2 7^ii/(|T-7yU,,J ^2 7^^ ^ W^(r:.,j,) (4.43) 

Since 6 = 6{f3, h,b) := Var^+ {a^;) > 0, from S439h . (lOTl) and (l4~43l) we see that 

IPp ((^1,/^) -mixing does not hold) ^ Pp^BS^,^, : ^ W^(r:«,j;) ^ 57^^^ 

^ W{Tx^y)^6j-''r''y (4.44) 

x: d{x) ^ L y&D^^i-^ 

For every x we can use the bound (I4.20D . so that (I4.44D yields 

?i, /i)-mixing does not hold) ^2 6-^ exp (-to5/'7"^i) . (4.45) 



Setting e.g. r = y^, = C log L with C sufficiently large we see that by the Borel Cantelli 
lemma we have (^i , r^^ )-mixing Pp-a.s. for all li large enough. This concludes the proof 
of the theorem. □ 

Remark. One may wonder whether the result of Theorem 14. 71 captures the true behav- 
ior of the logarithmic Sobolev constant in presence of the random realization of obstacles 
or whether, instead, it only provides a pessimistic bound. As we show below, as soon as (5 
is large enough (actually larger than the spin-glass critical point for the pure Ising model 
on tMH), in all the three cases described in the theorem, there exists a set of obstacles 
realizations of uniformly positive probability, such that for every uj ^ VLq the spectral gap 
and a fortiori the logarithmic Sobolev constant must shrink to zero at least as fast as L^'^ 
for some deterministic exponent C > 0. A quick sketch of the proof of this fact for the 
setting (a*) and e.g. 6 = 2 goes as follows. 
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Pick a vertex x G T{uj) with d{x) ^ ^ and denote by T^;^^ the finite sub-tree of rooted 
at X with i = 5{M + 1) log L levels, 5 « 1, M > 1 but (5M <C 1. Assume that T^^e, is free 
of obstacles, group together the sites of dT^/ into equal blocks according to their common 
ancestor in T^^i at level 6 log I/ + 1 and order the blocks from left to right. Then impose that 
all vertices inside the odd blocks are obstacles while all sites z inside even blocks are not 
obstacles and the corresponding Rz{uj) satisfies Rz(^) < £ (as usual e = e~^^). Because of 
Lemma 14.11 the obstacles realizations that obey the above specifications have probability 
larger than e~'^l^^'^^^^^'^l = e~^^" for a suitable constant c = c(p, j3) and a = 5{M+\) log 2. 
Since 5M <ti 1 the probability of finding a vertex x with the above properties converges to 
one as L ^ oo. 

Consider now the common ancestor y at level SlogL of two odd-even neighboring 
blocks. It follows immediately from the recursion (I4.2D and the assumptions we made on 
the even/odd blocks, that \Ry — 1\ ^ ^-cSM log l ^ suitable constant c. In turn, if M 
is large enough, that implies that the marginal of the Gibbs measure fj,L^^ on the finite 
sub-tree rooted at x with now 6 log L levels, has a bounded (independently of L) rela- 
tive density with respect to the Ising Gibbs measure on the same tree with free boundary 
conditions on its leaves. Since the latter has a spectral gap (and a fortiori a logarithmic 
Sobolev constant) smaller than 5^''(/5)'^i°g^ for some positive a(/3), we conclude that for 
the obstacles realizations satisfying the previous conditions, Csob(/^^^) ^ L^'^ for some 
C' = C'(A<5,6). 

5. Proof of Claims (1)-(4) 

The results of the previous section allow us to fill the gaps in the proof of Theorem ll.21 
We refer to section ISTD for the setting and the notation. 

Claim 3. Recall the definition of the random trees Tj 2^(1^) and the associated measure 
'^t'M' where z is a vertex at level i. We have to estimate min^ Cgap(z^^2T)- Since the 
b.c. above z can only affect this quantity by a constant factor (depending on /?), we may 
replace i^^'^^ by the measure iy^2e with free b.c. above z. At this point, for each z we are 
exactly in the setting of Theorem 14.71 with L = 21. As we have seen in the proof of that 
theorem (see (I4.45D ). we can prove that there exists C, < 00 such that 

Pp (csob(<2^) ^ ^ e-'' , (5.1) 

for all sufficiently large values of i. Since the number of z such that d{z) = £ is b^, the 
claim follows from the Borel-Cantelli lemma. □ 

Claim 4. Observe that it is sufficient to prove 

= +1) - ^r;,i(^x = +1) ^ (3&)~'') ^ . (5.2) 

Recall that i^^'^^ stands for the Gibbs measure on Tr^^2i{^) with — b.c. above Di{uj) and 
+ b.c. below D2i{uj). Observe that, by (14. 7D we have 

= +1) - ^r;i(^- = +1) = W^(rr-.,x) . (5.3) 

As in the proof of (I4.20D the expectation of the latter expression is estimated by 
'EpW{Vr^,x) ^ (2n)^/^, since x satisfies d{x,rx) ^ £/2 (see also Lemma l4~4l Lemma 1475] 
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and Lemma l4T6l) . Then by Markov's inequality 



1 



= +1) - c,2£(^- = +1) ^ 2 ^^^^ J ^ 2 " ^^-"^^ 

provided n is small enough. We turn to an estimate of the difference ^''^^^icrx = +1) — 

'^r^',oo(<7a; = +1). Let Ei{uj) denote the (Ising-model) event that there exists a path T in 
T{oj) joining the sets D^^^iuj) and 1)3^(0;) such that Uz = —1 for each z G F. If the sets 
Dbf^iuj) and D^i{oj) are not connected in T{ijj) we simply set Ei = 0. Observe that by 
monotonicity, for every uj we have 

^;:%c{<yx = +1 1 ED ^ v';:%,[ax = +1) (5.5) 

The reason for the above domination is that if there is no path connecting D5^{uj) and 
L)3£(a;) covered by — spins, then there must exist a cut-set of T{uj), fully contained be- 
tween level |£ and level 3£ covered by + spins, and conditioned on this event iyrjJ>o dom- 
inates i'~'to. We then have 

^ ^rSli''- = - ''rjU^e^rjU^x = +1 | E^ 

At the price of a /^-dependent factor we may replace i^rj^ with the measure ^ with 
free b.c. above the vertex r^. Now we are in the familiar setting of the previous sub- 
sections. To estimate ^oi^i)' suppose uj e is such that T{ui) contains a given path 

r = {xo,xi, . . . , Xfi} with d{xj^i) = d{xj) + 1, d{xo) = \i and d{xh) = 3i, h = i/2. Write 
{r = — } for the event {a^Q = a^^ = ■ ■ ■ = = — 1} and let qj{uj) denote the probability 
z^rt.ool-^ajj = -1 1 cr^j_i = -1). Clearly we have 

h h 

<,oo(r = -) = <,ooK, = -1) n ^ n • (5-7) 

Observe that 

qjiij) = e~^Rx^{u;)/{l + e-^Rx^{uj)). (5.8) 

As in the proof of Lemma 14.41 Lemma 14.51 and Lemma 14.61 we discriminate the vertices 
along r in good and bad vertices. We know that if z is good, then Rz{uj) ^ ue and 
therefore, by (15. 8t we have qj{uj) ^ u. As in (I4.24D we may then estimate 

<,oo(r = -)^«^"(") (5.9) 

where n{uj) stands for the number of bad vertices along T. Summing over all possible 
paths and estimating as in (I4.25t and (I4.33D we arrive at 

[i^tooiEi)] ^ b'' (2^)^ ^ \m-''e-' , (5.10) 
provided u is suitably small. From (15.61) . Markov's inequality yields 

(^x = +1) - '^rjU^x = +1) ^ ^(3^)"'') (5.11) 

This, together with 05. 4D . ends the proof of Claim 4. □ 
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Let us now turn to Claim 1 and Claim 2. Here the environment uj is given by a 
Bernoulli(p) configuration r] below level I and is deterministically free of obstacles up 
to and including level I, as prescribed by (I3.3t . 

Claim 1. Let Ei be the event that there exists a vertex x with d{x) = i, such that = — 1 
for every z G T^, i.e. if the root is connected to level £ by a path covered with — spins. As 
in (15. 5t and (15.61) we have 

^ - ^ . (5.12) 

Following fl5.7D and (15. 8D we estimate 

By monotonicity we have Rz{uj) ^ Rz{oj), where uj coincides with lj (and therefore with 
rj) below level £ and is given by a new (independent) Bernoulli(p) configuration i]' up to 
and including level i. We denote by Ep the expectation over the random environment uj. 
Here we can apply the machinery developed in Lemma IT4l Lemma |43] and Lemma 1431 
Namely, for a suitably small parameter u > 0, we can write e^^Rz{uj) ^ u for every good 
vertex z. Estimating as in ( 15. 9D and (15.11ft above we have 



x: d{x)=l 



Therefore 



^ {2uY ^ er^^ . (5.14) 



'^^,tm > e-^'j ^ e^'E, [fi+iE,)] ^ e^' . (5.15) 

Thanks to (15.121) and the Borel-Cantelli lemma, this implies the desired estimate. □ 

Claim 2. This is the same as the statement (I4.35D appearing in Theorem 14.71 with the 
difference that now the environment is deterministically free of obstacles up to and in- 
cluding level i. We can therefore repeat the argument used in the proof of Theorem 14.71 
and see that what has to be established here is a version of the exponential integrability 
(I4.20D for our new environment uj. The latter, in turn, relies on the bounds of Lemma l4.11 
Lemma 14.21 and Lemma 14.31 Since the ratios R are monotonic functions of the environ- 
ment, these estimates can only improve in the setting considered here and the proof of 
Claim 2 becomes a trivial modification of the proof of Theorem 14. 71 □ 

6. An extension to the hard-core lattice gas (independent sets) 

6.1. The hard core lattice gas. A configuration i] e := {0, 1}^'' is called an indepen- 
dent set if no two adjacent vertices are occupied, i.e. if rj^riy = for every couple x, y G 
such that d{x, y) = 1. We call the collection of all independent sets over the 6-ary tree 
T^. In the hard-core lattice gas model is the set of allowed configurations and each such 
configuration G O is weighted with the factor A'''' where |r/| stands for the cardinality of 
T], i.e. the number of occupied vertices in r/, and A > is the so-called activity parameter 
To define the Gibbs measure we use local specifications obtained by setting 

where ^4 is a finite subset of T^, T,r] e Cl are two allowed configurations such that = rj^ 
for all X ^ A, and = ^xeAVx- It is well known that the hard-core lattice gas model 
undergoes a phase transition at the critical activity Ac = — 1)^^^) (see e.g. 1127111011 ). 
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For A ^ Ac there is a unique phase regardless of the boundary condition on the leaves, 
while for A > Ac there are (at least) two distinct phases, corresponding to the odd and 
even boundary conditions respectively. The even boundary condition is obtained by 
occupying all the vertices at even depth from the root and letting all the rest unoccupied, 
i.e. 

g J 1 d{x) is even 
^ 1 d{x) is odd 

The odd boundary condition r° is the complement t° = 1 — r^. We use the notation 
/x^ = fj,Jf^ for the Gibbs measure on the tree of depth £ with even boundary condition. 
Similarly /x^ denotes the Gibbs measure with odd boundary conditions. We also write 
//^ = lim£_,oo fJ-i and fi° = lim^^oo When we need to emphasize the A-dependence we 
shall write /x^,/^^ in place of fi'^,fi°. Phase transition is reflected by the fact that, when 
A > Ac, the probability of occupation of the root differs for /^^ and 

6.2. The Glauber dynamics. The hard-core Glauber dynamics is the Markov process 
with Markov generator formally given by (11.11) with flip rates that are reversible w.r.t. 
the hard-core lattice gas Gibbs measure. As in the Ising model we restrict for simplicity to 
the heat-bath dynamics given by 



qx cr^ = l ^ ^ 

px a"" gCI, a^ = qx:= — ^ i P\ ■= j—j ■ (6-1) 

a"" 



Here cr G and represents the configuration a with the occupation number at x in- 
verted, i.e. {(T^)y = (Ty, for ally ^ X and {(t^)x = 1 — u^. In words, only transitions within 
D, are allowed and the transition a ^ occurs with rate px = X/{1 + X) if x is vacant 
and with rate qx = + X) if x is occupied. It is easily verified that detailed balance 
holds with this choice of rates. Moreover, for any finite subset A c T^, for any r G (7, the 
finite volume dynamics on A with boundary condition r is ergodic and reversible w.r.t. the 
Gibbs measure As for the Ising Glauber dynamics we can use the spectral gap and the 
logarithmic Sobolev constant to estimate the rate of convergence to the stationary distri- 
bution fi\. The corresponding definitions are exactly the same as in (12.61) . An important 
result of 111 911 is that the uniform bounds fl2.9D hold here if we replace with /i^, i.e. in 
the even phase one has exponential decay to equilibrium for all values of A. Of course, the 
same holds for the odd phase. 

6.3. Attractivity. It is essential for our approach that we can define a partial order on Cl 
such that the hard-core lattice gas and its Glauber dynamics become attractive. Let us 
write T'' as the disjoint union of even and odd vertices, Tcvon and Todd> where Tcvon := 
{x gT'' : d{x) is even} and Todd := {x G T^ : d{x) is odd}. We define the following order 
on Vl: 

I ^ X G Teven 

a ^T] <^=^ < (6.2) 

ycTx ^ Vx X e Todd 

A function / : ^ M is called monotone increasing (decreasing) if a ^ rj implies 
/(cr) ^ f{r]) if (a) ^ /(??))• We also write fi ^ u, for two measures on D,, when- 
ever fi{f) ^ u{f) for every monotone increasing function /. As in the Ising model it 
is straightforward to construct an order-preserving global path-wise coupling. Let erf '^'^ 
denote the hard-core Glauber process at time t, with start in the configuration ^ G il. 
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evolved in the region A with boundary condition r G 17. We may couple the pro- 
cesses {{crf^'^)t ^0, A c T^, G 17} such that the following relations hold: for any 

yl C S C T^ any ^ -< r -< r' and i ^ 

(T^' ' -< fj^ ' ' (6.3) 
fjj ^ fjj ' ^ fj^' (6.4) 

These relations also imply the following monotonicity properties of the Gibbs measures 
and the associated FKG-property (see (I2.1I) - (I2.2I )): 

(i) for any A cT'' the map ij fJ^'j^if) increasing; (6.5) 

(ii) ^"^Q ^ jjJ^ whenever Ac B. (6.6) 

6.4. Results. Replacing /i+ with we may define the sets IIq,, a G (0, 1), just as in 
Definition ll.il It is not difficult to check that Lemma IZD and therefore Corollarv l2 .21 hold 
in the present setting as well as in the Ising case. The same applies to Corollav 12.31 and 
Lemma IT4I 

We need to introduce the hard-core analog of the Bernoulli measures Pp. We call Vp^x, 
p G (0, 1),A ^ the probability measure on obtained as follows: we first assign oc- 
cupation numbers on Tcvcn according to the Bernoulli(p) probability Pp. This gives a 
configuration r/ on Tcvcn- To obtain a legal configuration (in Vl) we may now occupy only 
those vertices in Todd that are at least at distance 3 from 77. Call c Todd this set of 
available vertices. Finally put rf^ = 1 with probability p\ = A/(l + A) independently for 
every x e A^j. 

A simple coupling argument shows that, for every A > 0, fp a ^ /"a soon as p ^ pa- 
In particular, the argument of Lemma l2.4l shows that I'p^xi^a) = l, for every a > 0, for all 
p ^ Px. Our main result for the hard-core lattice gas is stated as follows. 

Theorem 6.1. 

a) For every b ^ 2, there exists p < 1 such that for all A G (0, 00) we have ^{Qa) = h for 
some a = a{X, b) > 0, for any initial distribution v such that v ^ Vp^\ . 

b) For every p> \, there exist bo e'N and Aq G (0, 00) such that for b ^ bo, X ]^ Xq we have 
i^(Oa) = 1, for some a = a{X, b,p) > 0, for any initial distribution v such that u ^ Up x ■ 

6.5. Sketch of proof of Theorem 16.11 Theorem 16.11 will be proved with the same ar- 
guments used in the proof of Theorem 11.21 Below we point out the necessary (rather 
obvious) modifications. 

The first observation is that in view of the monotonicity of Ua, the domination Up^x ^ /^a 
for p ^ px, Lemma |Z41 allows to replace the statements in the theorem by 

a*) For every b ^ 2, there exist p < I and Aq < 00 such that for all A ^ Aq we have 

i^p^xi^a) = h for some a = a{X, b) > 0. 
b") For every p > ^ there exist bo and Aq such that for 6 ^ 60, A ^ Aq we have Up^xi^a) = 1 

for some a = a{X,b,p) > 0. 

To repeat the argument of section|3lwe need to introduce the notion of the environment 
of obstacles. A realization of the environment is described by a; G 17 with the following 
interpretation: x G Toven is called an obstacle if u)-^ = and is said to he free if lVx = 1. 
Similarly x G Todd is an obstacle if oj^ = 1 and is free if ujx = 0. Note that x is an 
obstacle in oj iff = t°. As in the Ising case lv determines the tree T{uj), i.e. the largest 
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connected component of free vertices containing the root. We write B^o for the set of r e 
such that Tx = for every x ^ T{uj). The hard-core model in a given environment uj is 
then obtained as before: for every yl c we write fi^ ^ for the measure fJ-^^j^^^^y where 
T ^ Bu). The same reasoning apphes to the dynamics and we may use, as before, pt,uj{0 for 
the expected value at the root of the occupation variable under the dynamics a^^ among 
obstacles with starting configuration ^ e 
We then observe, as in (13. 2D that 

where pt{uj) denotes expectation at the root w.r.t. the dynamics in infinite volume without 
obstacles with starting configuration uj. We then define the environment uj = uj{ri,i) as 
in (I3.3D . where of course the + configuration is replaced by r^. We now proceed exactly 
as in (I3.5D . Moreover, we may repeat the estimates of the three terms there without 
modifications. What is crucial is that the technical estimates isolated in Claims 1 to 4 can 
be established for the new setting. A discussion of the point is given in the next subsection. 

6.6. Technical estimates. To prove the Claims 1 to 4 for the hard-core model one needs 
to adapt to the present setting the analysis developed in section|4l One defines the ratios R 
and the associated weights in a similar way here, but the recursive relations involved in 
the proofs of the main estimates are model-specific and require a separate investigation. 
We will not provide all the details here since there is no truly new ingredient. However 
we give a sketch of the basic computations on the ratios R to help the interested reader in 
reconstructing the needed claims. 



Estimates on R. Let R be defined by 



f^ZicTr = 0) 



A simple computation gives that if = 1 we have 



where Xi, i = 1, ... ,6 denote the children of the root and Rx,{^^) is the corresponding 
ratio, given as usual by the rule Rx{i^) = R{6xUj) [9xUJ being the environment shifted by 
X G T^) . The crucial estimate (I4.10t is now replaced by 

R^Vx) ^ 5 . (6.9) 



Here a denotes the probability i/p a conditioned to have uJr = 1 and 5 is a small param- 
eter to be fixed at a later stage. The following bound is the analogue of Lemma RTTl in the 
present setting. 



Lemma 6.2. For any 6 > 0, b ^ 2, there exist pq < 1 and Aq < oo such that (f6.9P holds for 
all p ^ po and A ^ Aq. 

Proof. We shall give the proof only in the case 6 = 2. For any integer £ we may define the 
ratios R^{uj) w.r.t. as in (14. IIP . We set 



JV,A (i?^ ^ Va) (6.10) 
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Let xi,X2 denote the children of the root and call yi,y2 and 2/3,2/4 the children of xi and 
X2, respectively. Observe that the event E that ujx, = 0, i = 1, 2 and ujy. = 1, i = 1, . . . ,4 
has z/p A probability at least 1 — 4(1 — p) (since it suffices to occupy all y/s to automatically 
free the Xj's). Moreover for a; G £^ we have 



R^u) = A 1 + zrr^-— ,rr-r^ U + 



A \-V. A 



[1 + Rl^{co)){l + R^^ico)) J V {I + Ri,{u;))il + R^^iu;)) , 

(6.11) 

Note that Ry^iu)) = R^''^{9y^Lo). Suppose that Ry^iuj) > \/A, i = 1,2,3. Then the above 

formula shows that for A sufficiently large, the condition R^{uj) ^ a/A forces Ry^{i^) ^ 3. 
Reasoning as in (I4.14D we see that 

^ 4(1 -p) + 6qj_2 + 4i7p,A (< ^ 3) . (6.12) 

Using again (16.111) we see that if there is only one yj with R^. {uj) ^ a/A it is impossible to 
have R^{uj) ^3. It follows that 

^P,x (4i ^ 3) ^ 4(1 -p)+ 6(/|_4 (6.13) 
Putting these estimates together and using qe^^ ^ g^_2 we see that 

g£ ^ 12(1 - p) + 30g|_2 . (6.14) 

The conclusion now follows from (I6.14D just as in the case of (I4.15t because of the even 
boundary condition. □ 

7. Open problems 

We conclude by discussing an interesting open problem. Back to the Ising case with 
h = 0, let us take as initial distribution for the Glauber dynamics the symmetric product 
measure ¥1/2 that for shortness we denote by u. 

A first non trivial question is whether the law of the Glauber dynamics i/Pt converges 
to a Gibbs measure as t ^ 00. In Z'^ it is well known that this is the case (see e.g. 
II16II ) because uPt is translation invariant; unfortunately the Lyapunov function techniques 
behind the proof do not seem to apply on the tree because of the large boundary/volume 
ratio. 

In the uniqueness region /? ^ /?o it is not difficult to check that vPt converges weakly to 
the unique Gibbs measure as t ^ 00. More interesting is the interval /? G {Po,Pi)> where 
Pi is the spin-glass transition point discussed in section 11.11 Here the situation is more 
complicate due to the presence of infinitely many extremal Gibbs states. 

If we recall our first characterization of Pi, it is not unreasonable to conjecture that uPt 
will converge to the (extremal) free Gibbs measure In fact, if we imagine that the 

single site Glauber dynamics is replaced by a block heat bath dynamics as in [ 1 ] then, at 
least for small times, each update of a block (say a large but finite subtree) replaces the 
Bernoulli product measure u inside the block with a finite Gibbs measure close to fi^'^^^. The 
case P ^ Pi should be even more complex and one can conceive that the dynamics and 
coarsening of clusters of spins with opposite sign, present in the starting configuration, 
will play a significant role as in the P = +00 case 113 111 . 

Although we have no clear answers to any of the above questions, we do have some pre- 
liminary "concentration of measures" results that bring some support to the conjectured 
behavior in the intermediate regime (/?o. Pi)- 
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First we show that for any local function / 



^(^?; \Pt{f)iv) - i^Ptif)\ > e 



ct \ ^ g-ce 



ct 



for some c > 0. 

In other words not too large fluctuations in the starting configuration rj are completely 
washed out by the dynamics. In particular, for any local function / which is odd w.r.t. a 
global spin flip, 

limPt(/)(r?)=/-^(/) = z^-a.a. r/. 

t^oo 

Secondly we derive a stability result that can be roughly formulated as follows. Let i> be 
a perturbation of u such that the relative entropy between u and i> restricted to the first i 
levels does not grow faster than (b^Y, 6 '^1. Then for any local function / 

lim \i)Pt{f)-iyPt{f)\=0. 

t—fOO 

We now formalize what we just said. 

Proposition 7.1. For h = and /3 < /3i there exists a positive constant c > such that, for 
any function f depending only on finitely many spins and any t ^ : 



^[m \Ptfiv) - ^Ptf\ > e-^*j ^ e-'f' . (7.1) 
for a suitable constant c/ > depending on f. 

Proof We are going to use standard Gaussian concentration bounds 111 511 for the measure 
z/ of the form: 

i/(r/; |F(r?)| r) ^ 6"^ (7.2) 
for any mean zero function F with unitary Lipshitz norm 

\\F\\h^ := Yl mr^i-pml- 

Therefore (17.11) follows if we can prove that for some a = a{P) > 

||PJ||Lip<C/e-"*. (7.3) 
for a suitable constant C/ > depending on /. The basic tool is coupling along the lines 
introduced in |[T||. Recall that tanh(/?) < l/\/6 for any /? < Thus we can always choose 
A G (tanh(/3), (6tanh(/3))~^) in such a way that < 1. Given two configurations r/,^ 
that differ in finitely many points, define their weighted Hamming distance as 

X 

Then a key result of Q combined with an unpublished paper of Peres and Winkler (see 
section 4 of Q) shows that under the natural coupling of the Glauber dynamics started 
at r] and ^ 

E(dA(a?,4)) <Ce-^*dA(r?,e) (7.5) 
for suitable positive constants C, c. Therefore 

y 

^ A-^(^) Wfiv^) - /(r?)||oo)e-^*A'^(^) = Cf e-^X"^^^ (7.6) 
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Since bX"^ < 1, the sum over x of the square of the r.h.s. of (I7.6D converges and (I7.3D 
follows. □ 

Corollary 7.2. 7n the same setting as above, let v he a prohahility measure on fl and let 
ui, vi he the marginals on $7^^ ofu and u respectively. Then there exists 6 = 6{P) such that, 
if Entj.(§^) < b^^ for all I, then 

\iui\DPt{f)-uPt{f)\=Q. 0.7) 
Moreover the limit is attained exponentially fast. 

Proof Let k = k{(3) be so large that, with £ = kt, for any large enough t 

\\Pt{f)-e"''r\f)\\^^e-\ 

where C^l^'^ stands for the generator of the Glauber dynamics in with free boundary 
conditions. Standard results on finite speed of information propagation show that such a 
k exists (see e.g. [18]). Let now c/, c be the constants appearing in Proposition l7.1l and let 

At be the set of configuration {q; \Pt{f){'n) - vPt{f)\ ^ e^"^*}. Then, by setting := 

\m{f) - i^PtU)\ ^ 2e-* + \D,e"''r\f) _ u,e"''r\f)\ 

= 2e-* + |z.([/i,-l]e<°^(/))| 

^ 4e-* + 2e-^* + e"^/^" + \\fUv{hdA,) (7.8) 

It remains to bound z/(/i^lAt) and this is easily accomplished using the entropy inequality 
together with Proposition 17 11 and our assumption on Ent;^(/i£). For any A > 

u{haA,) ^ ilog(z^(e^i-4*)) +\^^tu{hi) (7.9) 

^ - log(l + (e^ - 1)6-'=/^") + — (7.10) 

If we now choose A = \cfe^ and 5 < -j^^^ we see that u{hilAt) tends to zero as t ^ oo 
exponentially fast. □ 
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